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I.  Olkin,  M.  Sobel  and  Y.  L.  Tong 

Stanford  University;  University  of  California, 
Santa  Barbara;  University  of  Nebraska 


Abstract 

The  usual  formulation  of  the  problem  of  selecting  the  best  of 


k populations  has  been  studied  from  the  point  of  view  of  designing 
an  experiment  that  will  guarantee  a certain  probability  of  making  a 
correct  selection  (PCS).  However,  in  retrospective  studies  the  sample 
has  already  been  taken  and  we  wish  to  estimate  the  probability  of 
making  a correct  selection  using  this  data.  The  present  study  deals 
with  obtaining  a class  of  upper  and  lower  bounds  for  the  true  PCS  and 
its  estimator  using  the  methods  of.  majorization.  We  also  study  the 
large  sample  behavior  of  the  estimator  of  the  PCS,  which  turns  out  to 
be  asymptotically  normal. 


Int roduction. 


Consider  an  independent  observation,  X.  from  k populations 

1 -3 

F(x,0. ),  i = 1,  ..,k,  3 = l,...,n.  The  experimenter  wishes  to  select 
the  population  associated  with  the  largest  parameter  0^]*  To  do  this 
we  define  an  appropriate  statistic 


Y.  = Y (X. . , . . . ,X.  ) , 
i il’  ’ m ‘ 


and  use  the  decision  rule  "designate  the  population  corresponding  to 
the  largest  Y.  value  as  the  one  with  the  largest  0-value."  If  the 
c.d.f.,  G(y,0^),  of  Y^  has  the  form  G(y-6i),  then  the  probability 

of  a correct  selection  (PCS)  is 


(1.1) 


r kri  1 

TT  Gty+b.^  dG(y)  , 


where  P ■ < •••  < 0f.  , „ 5.  = 0-  - 0r . , • 

LU  - - [kj  i [kj  [ i ] 

In  the  usual  methods  employing  the  "indifference  zone " approach, 
as  introduced  into  ranking  and  selection  by  Bechhofer  (195*+ )>  we 
determine  the  smallest  fixed  sample  size  n needed  in  order  to  satisfy 
a preassigned  probability  requirement.  This  led  to  a consideration  of 


least  favorable  configuration,  which  frequently  takes  the  form 

s 6n  and  suggests  an  inequality  for  Of ), 
where  we  "move  downhill"  in  the  sense  of  majorization  to  the  case  of  equal 
&. -values,  regardless  of  whether  the  common  value  is  specified  or 

not.  Because  this  integral  is  often  difficult  to  evaluate  (especially 
in  high  dimensions),  we  are  interested  in  both  upper  and  lower  bounds 


•» 


that,  hold  for  any  true  configuration,  without  regard  to  any  least 
favorable  configuration. 

In  the  present  paper  we  consider  the  problem  of  estimating  the 
true  PCS  by  using  information  provided  by  the  sample.  This  is 
extremely  important  for  those  problems  in  which  we  are  carrying  out 
a retrospective  study  — rather  than  designing  an  experiment. 

If  Y is  a consistent  estimator  of  6,  then  we  may  let 

Yr - i < ...  < Yr,  i be  the  ordered  Y-values,  and  define 

UJ  - - [kj 


&i  " Y[k]  ' Y[il  ' 


1=1,.  . ,k-l  . 


Then  6 = will  be  a consistent  estimator  of  5,  and 

. A A 

, . . . ),  will  be  a consistent  estimator  of  a(&  , ...,5  ). 

M A™  J-  " -L  K-_L 

As  a consequence,  we  wish  to  further  study  the  behavior  of  the 
integral  . /}  or  ^,s  as  a function  of  the 

configuration  of  the  parameters.  In  particular,  we  are  interested  in 
the  following  question:  how  does  the  PCS  behave  if  some  of  the  5-values 
are  replaced  by  their  average:  This  problem  is  studied  in  the  next 
section. 


Pemark:  The  above  development. was  based  on  a location  parameter  family. 
However,  a similar  analysis  holds  for  scale  parameter  families.  In  this 
case  the  true  PCS  (l.l)  is  replaced  by 


(1.2) 


a 


(&!>..., &k_l)  = J ^ j J G(&iy)j  dG(y)  , 


2 


* 


*1* 

r 


i 


where  in  the  case  of  the  largest  scale  parameter  we  define  6.^  = ®[k]/e[i] 
(i  = 1,2, . . . ,k-l ) and  hence  6^  > • • • > 6 2 > 1.  For  the  case  of  the 
smallest  scale  parameter  we  define  ^ then  (1.2) 

and  the  above  inequalities  on  6^  still  hold. 


2.  The  Main  Bounds. 

The  key  motivating  point  is  that  under  certain  circumstances  we 
can  approximate  ®{b)  by 


® • • • f t’p  ) = ® • • • > f 2_, 


where  Z p^  = p = k-1  and  the  t^'s  are  certain  averages  of  the  b's. 
If  the  approximation  sign  is  an  inequality  sign,  then  we  are  "moving 
downhill"  by  averaging  and  further  averaging  will  move  us  more  downhill. 
This  suggests  the  use  of  majorization,  and  we  show  how  this  arises. 

We  now  develop  the  main  result.  For  any  fixed  r < p,  let 
a1,...,aj.  be  distinct  numbers  satisfying 

0 = aQ  < < a;  < • • • < s p , 

p = a -a  v = l,...,r,  and  define  the  partition 

v v v-1 


&(A)  = (£^,...,6  , ^a.+l’  ’ * * ’ Sa.  1+l,"*,Jp' 


r-1 
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We  also  define  the  means  within  each  cluster: 


and  generate  a vector  of  averaged"  o- values : 


where  A denotes  the  set  A = (a^,,..,a  } . 

Let  A c B mean  that  B is  a refinement  of  A-  A question  that 
arises  immediately  is  whether  we  can  compare  the  O -values  of  two 
partitions  A and  B if  A c B? 

Theorem  1,  If  leg  G(y)  is  concave,  and  A c B,  then 

(2.1)  a(61,...,&p)  <a(61(B),...,&p(B))  ^a^CA),...,^  (A))  . 

The  proof  of  this  result  is  given  in  the  next  section.  We  note 
that  the  cr's  in  (2.1),  though  written  as  a single  integral,  are 
actually  p-fold  integrals.  However,  we  would  also  like  to  have 
t * bounds  that  reduce  the  dimensionality.  To  accomplish  this  we  consider 

i> 

t the  following  integral.  For  any  partition  A,  define 


*✓  *. 
► • 

. ’it 
* 


«> 

c 

V *. 

v*  * . 


i 


(2.2)  P(B1(A),...,B  (A)) 


a_ 

oo  2 


-oo  a^+1 


G (y+&i  ) dG(y) 


-[£  n°(yt6i)dG(y)]  [/. 

...  f f h G(y+B.  ) dG(y) 

J-  oo  a '+i 


Here  the  reduction  in  dimensionality  is  in  terms  of  the  size  of  the 

intervals  a - a . (v  = l,2,...,r).  Note  that  when  r = p,  then 
v v-l 

P (B^j  • • • > 6p ) = Gt  (6^/ « * « ) * 

Theorem  2.  If  Ac  B.  then 



(2.3)  P(B1(B),...,B  (B))  < P(61(A),...,B  (A))  < a (^ . . . ,&p)  . 

Of  particular  interest  are  the  extremes  obtained  by  the  partition 
A = (1, ...,p},  in  which  case  we  obtain 


(2.4) 


P r00  r 00 

|T  j G(y+B1)  dG(y)  <0(61,...,B  ) = / f|  G(y+6. 

1 J -00  * J -00  1 


) dG(y) 


Remarks : (i ) We  should  note  that  Theorem  2 is  due  to  Kimball  (1951 )• 

(ii ) If  G is  a normal  distribution,  then  the  right-hand  side  of 
(2.3)  can  be  regarded  as  the  integral  of  a multivariate  normal  distribu- 
tion with  variances  1 and  correlations  p = l/2.  This  integral  has 
been  considered  by  Milton  (1970,  1972).  If  G is  normal  then  the  left- 

hand  side  of  (2.4)  is  easily  shown  to  be  -fr  G(B«//2),  so  that 
integration  is  required. 


no 


i=l 


? 


5*  ' v 


(iii ) If  G(y,0)  = G(y /q),  with  G(0)  = 0,  is  a scale  parameter 
family,  then  we  define  geometric  means  instead  of  arithmetic  means. 

For  any  partition  A we  have 


/ al  \1/al 

» 1 S2  ) 

\Y  6i)  ' 

HttO 

l/ (a2“ai ' 


6 = 

r 


t h 


r-1 


L/ 


The  results  of  Theorems  1 and  2 hold  with  the  obvious  changes. 


Proof  of  Inequalities  via  Ma.jorization. 

We  first  review  several  facts  concerning  the  partial  ordering 


majoriza^on.  A vector  u^, . . . ,u  is  said  to  majorize  a vector 


^,...,v  , written  u V v,  if,  upon  a rearrangement  with  u^  > 


> v 


m 


m 


^ £ J> 

'■<  - ’ * ' > we  have  Z u-i  > X vs  > m = 1, . , . ,p-l,  ^ u = v . 

x p 1 1 “ I 1 l1! 


The 


class  of  (differentiable)  order  preserving  functions  f(x^,...,x  ) are 


called  Schur  functions,  or  Schur -convex  functions,  and  are  defined  by 
the  property 


1 J ’ 


for  all  i, j = 1, . . . ,p 


(For  the  Schur-concave  functions  we  merely  turn  the  inequality  around. ) 
Then  one  characterization  of  the  Schur  functions  is  through  the  fact  that 
u ) v if  and  only  if  f(u)  > f(v)  for  all  Schur-convex  functions. 


Finally,,  there  is  another  characterization  in  terms  of  doubly- 

stochastic  matrices  P = (p.  .),  p.  ..  > 0,  L p.  . = £ p.  . ~ 1>  namely, 

J 

u / v if  and  only  if  v = uP  for  some  doubly  stochastic  matrix  P. 

This  says  that  one  "moves  downhill'  from  u to  v (in  the  sense  of 
majorizaiion ) by  taking  a convex  combination  of  t ments  of  u. 

We  need  the  following  fact- 

Fact.  [Marshall  and  OLkin  (197*0].  If  log  G(y)  is  concave,  then 

p r oo  ,p 

||  G(y+£. ) is  S cnur -concave.  Consequently,  a(&)  = I j [ G(y+E^)  dG(y) 

1 J - 00  1 

is  S chur-concave. 

The  proof  of  Theorem  1 is  a direct  consequence  of  the  above 
facts.  Recall  that  in  (2.1)  we  have  that,  if  Ac  B then 

a(d ) < a (8(b))  < a(d(A))  . 

We  need  orny  prove  the  left-hand  inequality  since  the  right-hand  inequality 
is  obtained  by  a repeat  ed  argument  where  cu&fB))  is  used  in  place  of 
O(b)  and  ol  [b  A))  is  used  in  place  of  oc(6( B)). 

It  is  immediate  that  6 1-  &(B)  since  &iB)  -•  &P,  where 
P block  disg ! P. , . . . ,P  ) is  a doubly-stoehast.ic  matrices  with  blocks 
P on  the  diagonal  consisting  of  l/^a  ^-a^  1 ).  This  majorization 
pi ^s  the  Sonar  concavity  proves  Theorem  1, 

Theorem  2 follows  from  an  inequality  of  Kimball  [1951]  and  to 
apply  this  it  is  sufficient  that  the  distribution  functions  G(y+&^) 


all  be  positive  and  monotone  in  the  same  direction  for  all  real  y. 


uraximabiuns . 


We  now  make  use  of  the  bounds  to  approximate  O'-(E  ).  A3. though 

a (f, ) can  be  computed  exactly  with  the  help  of  quadrature  methods,  the 
necessary  programming  is  not  always  available  and  the  computations  may 
be  costly.  Thus,  it,  is  desirable  to  be  able  to  approximate  the  integrals. 

As  an  i. lustration  of  the  bounds,  consider  a numerical  example  for 


the  case  of  normality  with  common  cr  - 1,  and  the  location  parameter 
model  with  & - (3.2,  2.7,  1.5-  1.9,  1.7);  we  wish  to  evaluate  a(b). 

The  exact,  vauue  is  Ct  = 0.80l6.  To  compute  our  smallest  lower  bound 

5 

note  that  / G(y+&)  dG  [y  . so  that  we  obtain  j $f& . j/2)  = O.Jbh^ 

- 1 1 

The  largest  upper  bound  is  / ^(y+b)  d$(y)  = 0.8483,  obtained  from 
Milton  (1970). 

Wre  now  consider  some  refinements.  Suppose  we  use  A = {l,  3,  5) 
so  t.hw  we  cluster  into  3 groups  {3 .2),  {2.7,  2.5},  {1.9,  1*7),  then 
the  lower  bound  & (&  (A) ) in  ) is  given  by 


f $ (y+3.2 ) dt>(y)  • <J> (y+2.7 ) ®(y+2.5)  d®(y; 

* / $ (y+1.9)  <t(y+1.7)  d«t(y)  , 


vr.ich  may  be  obtained  from  the  Milieu  (l.TO)  tables  to  be  0.7698. 
But  if  we  now  approximate  (U.l)  by 


(>.2) 


I $ (y43. 2 ) d<J>(y)  / (t*1*  (y<-2.6 ) d<f(y)  / 0 (y+1.8)  dl'(y)  , 


then  we  obtain  the  result  0.771?. 


8 


f 


If  we  use  fewer  clusters,  namely,  A = (3,5),  yielding  (3"2,  2.7,  2.5) 
(1.9,  1.7)  then  we  obtain  the  approximation 

(4.3)  / ^(y-t-2.8)  d^(y)  / e2(y+1.8)  d$(y) 

to  yield  a value  of  0.7802. 

It  is  intuitively  clear  that  if  we  can  take  advantage  of  the  close 

/N 

proximity  of  any  two  5 values,  then  the  worst  case  for  getting  close 
bounds  and  approximations  may  be  the  equally  spaced  case.  We  have, 
therefore,  attached  at  the  end  of  this  paper  a table  (Table  l)  for  the 
normal  distribution  under  equally  spaced  configurations  with  a common 

spacing  d.  The  table  gives  the  numerical  values  of  the  integral 

P _ 

/ TT  <J>(y+id)  d<f(y)  = P[Zi  < id/,/2,  i = l,.,.,p]  for 

d = 0.00(0.02)0.20(0.10)4.00  and  k = 2(1)10,  12,  15  and  20; 
where  <J>  is  the  standard  normal  cdf. 


5.  Limiting  Distribution  of  the  Estimator  of  the  PCS  for  the  Normal 
Location  Parameter  Model. 

In  this  section  we  obtain  the  asymptotic  distribution  of  the 


estimator,  a(&),  of  the  true  PCS  a(b)  when  the  underlying  distri- 
bution is  normal.  Although  results  for  the  location  parameter  case 
are  obtained,  similar  results  may  be  obtained  for  scale  parameter 
families. 

Let  Y = (Y1,...,Yk'  and  Y = (Y[]_  j > • • • >Y[kj  where 
Y[l'  — < Y[k'  are  the  ordered  Y-values.  Similarly,  we  define 


9 


r 


IU 

C 


<- 


(01..  ...,6k)  and  6=  (0[1]>‘«*>erk])1  Finally  we  let 


= (0  : ® j i | < ” ' < 6[k , J with  strict  inequalities  between 


d-vaiues.  For  any  two  vectors  u and  v we  write  u < v to 


mean  u.  < v.  for  all  i. 
1—1 


New  suppose  Y^-6^))  - --'11(0,1).  We  first  obtain  the 


asymptotic  distribution  of  Y,  then  of  estimators,  , of  £, 


and  finally  of  «(&).  In  the  following  development  <r  is  assumed 
to  be  known,  and  we  take  it  to  be  unity. 


Theorem  3.  As  n ->oo,  tQri  (Y-e))  — »*T\,(0,I  ) for  every  6 e 0, 


where  I is  the  identity  matrix  of  order  k. 
K 


Proof.  The  result  follows  by  bounding  the  right-hand  side  of 


(5.1) 


P^(  -/n  (Y-9)  < z)  - TT  $(z^)| 


' |p  {-/n  (Y-e)  < 2,  Y = Y)  - P_  {-/n  (Y -C'  < zl 

e ~ e 


+ P fyn  (Y-e)  < z,  Y £ Y)  + |p^  (y'n  (Y-?)  < z) 

e e 


K 

TT  : 

1 i 


for  every  z = (z. ,...,z  ) and  every  e e n . 


Let  > 0,  e0  > 0 be  arbitrary  fixed  real  numbers  such  that 


< min<'er.  - - 6r.  i i )•  Since  Y is  a consistent  estimator  of 
1 - i ' L l ■ 1 1 -1 1 


and  the  YTs  are  independent,  it  follows  that  for  every  P c-  *2 


and  every  z - (z^,...,z,  ) there  exist  numbers  N1  s (c^,  l0  ' 


and 


N-  s N,.  (e  ,,  z ) such  that  for  n > N. , P(Y=Y)>P{|y.  -6.  ]<€.., 

1 6 ~„e  1 1 k 1 

l,...,k}  > 1-tp,  and  for  n > N_t,  !P  (y'n  (Y-~  < z)  - 7 T <f(z.  )!< 

- d d £ - 1 1 - e! 
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Here  is  small  enough  so  that  Iy^-P^i  — €1  implies  that  Y = Y, 

so  that  the  first  difference  on  the  right  side  of  (5.1)  is  zero.  Then 
for  n > maxtN-^rt,},  (5-1)  is  bounded  by  2e0,  which  completes  the  proof. 

As  a consequence  of  Theorem  3,  we  have  the 

Corollary  h.  For  every  6‘  e Q,  £ (J~n  (t  -l  ))  is  asymptotically  normal 
with  mean  vector  zero,  variances  2 and  covariances  1. 

Now  consider  the  PCS  for  the  normal  case 


/“  J0_ 

TT  ®(t  +v/n  Bp)  dfl>(t)  , 

■00  1 


and  its  estimator. 


~ ~ ^ r00  P _ ^ 

a(B)  = a (b2 ,...,£>)  = |[  ®(t  + -/n  ^ 

J —00  1 


) d4>(t ) . 


Theorem  t.  £(Jn  [a(&)-a(&)])  H(0,vn),  where  vn(&)  = ^ d2  + (f  d^2, 


/-  /■/  n 5. 

(5-2)  d.  = ^ qp  — ± 
•/2  -/2 


r)£  J[ 


Proof-  A Taylor  expansion  about  6 = e>0  = (6^,...,?^)  yields 


^ (a(6)-a(60))  -ei0)  2 fe^io)2 

1 » x Af  . 


I &=SC 


a2a(f) 


+ 7 (£  .*  )(P  -f  ) A..,/ 

2&  1 10  * 


s wx(n)  + w2(n)  , 


where  6.  < b.  < b. . Corollary  J*  and  the  boundedness  condition 
1 — iO  — 1 


j | < c/6^  implies  that  £ (iv  (n 


b ( f v / . Consequently, 


the  proof  will  be  complete  if  we  show  that  plim  w.  (r:  - . To  see 

c 

2 

this,  we  use  the  fact  that  x<p(x/  and  x q>(x ' are  bounded,  s:  that 
there  exist  positive  constants  c^,  c^,  c..  such  that 


Sga(&)[  . C1  1 5Pa(&)|  ^ C2  , c3 


*1  ' fi 


lot,  - 5. 
i j1  it, 


If  we  define 


Kn)  .4j  £ («1-.i0)2  X .■/§  J 


c2  c, 

~7 ' ~ ' 


then  plim  U(n)  - 0.  Consequently,  if  Q = > £*,  i = l,.,.,p}  , 

for  any  b*  < b^}  it  follows  that 


(5-3)  P{  |w2(n)|  > e)  = p(|w2(n)!  > e,  C.)  + P{|w?(n,'|  > e,  G) 


< p{un  > e)  + p(a) 


But  for  fixed  e > 0 and  sufficiently  large  n,  each  term  on  the  right- 


hand  side  of  (5.3)  is  < e. 


The  asymptotic  variance  vn  defined  in  Theorem  1 is  quite 
complicated.  We  may  use  & in  v (6)  to  yield  a confidence  interval 
for  the  PCS,  a(&),  namely, 


12 


/n  (a  (8)  - a(S)| 


n(o.i)  • 


vn(6) 


However.,  we  may  also  wish  to  approximate  the  asymptotic  variance  v (8) 
or  its  estimate  v (&).  We  again  .use  majorization  to  obtain  some  bounds. 
To  do  this  we  require  the  following  facts 

Fact  1.  For  fixed  i,  and  any  z,  the  function 


TJ  ®|  -p  + /n  6 
dVi  1/2  J 


7n  6. 

T l 


is  a Schur-concave  function  in  , , &.  ). 

1 i-l'  l+l  p 

This  follows  from  the  result  of  Section  3. 

Fact  2.  For  fixed  i and  8,  > • * • >8  , 

1 - - tr 


~ &i  &.  6. 

6 = r'l  " ~ Si~l  " ~2  ' &i+l  ‘ ~ 


B. 

- £p 


^ (p^T  & “ 2^17  si'---'  5^1  6 


2(p-l)  i ~ L 


BJ  * BT  , 


- £ , 

where  B = 2*  /P  • 

1 

To  prove  this  result  we  note  that  the  elements  of  6 are  means  of 

L 

the  elements  in  8 . 


Fact  3*  For  fixed  i and  &„>•••>&  , 

— 1 - - v 
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B.  B B. 

6 s|&i  fi-i  - - * &i+i 


• e. 


, b B 

< |P&  “ (P“2)&p  “ 2 &p  " T - 4 


To  prove  this  majorization  we  first  note  that  p&  - (p-2  t - — £ 

e * p 2 : 

> & - &p/2.  Further,  the  sum  of  ail  the  elements  in  each  vector  is  equal 
to  p&  - 4(p+l)&^*  Consequently,  we  need  show  that  for  m l,...,p-l. 


B.  p mB . 

“(6p  - “?)  < l V~T 

y p-m+1  J 

jA 

But  this  inequality  is  equivalent  to 


p-m+j 

j#. 


6j/(m-l)  > B - (Bp-Bj)^  , 


which  clearly  holds. 

As  a consequence  of  the  majorization 


&L'\  bU  , 


we  have  that  for  every  z and  every  i 


i?"1  -1  + P&  . (P^i;  .1  s D (t n 

P-1  2Tp-i)  i|  2lt' 


|>lt,j^Z+t/np&  - (p-2)  vn"  &p  - &i)  ^ | ~ + Vn  &p rr--J  Dj 


It. 


lU 


Remark.  The  upper  bound  is  a monotone  decreasing  function  of  b^, 


and 


the  lower  bound  is  a monotone  increasing  function  of  b 


Consequently, 


if  5^  > b *,  we  may  obtain  the  bounds  where  b*  is  inserted  in  place 

of  b^  in  D2(t)  and  in  place  of  b^  in  D^Ct).  Call  these  functions 

D.(t),  j = 1,2.  From  the  definition  of  v , we  see  that 
j n 


C 

» 

e,  v 

'r  ' . 


EJ-  H(B)  [/  D1(t)  d$(t)]2  < vn  < ^ H(B)  [J  D2(t)  d®(t)T 


where 


H(8)  = ^ exp(-n£2/2)  + | £ exp (-nb2/U )J  . 

The  main  point  is  that  under  certain  circumstances  we  may  approximate 
vn  by  n K(b)/hn,  which  is  small  for  large  n.  This  holds  when  the 
dispersion  of  the  b-values  is  not  too  large,  so  that  b,  < 2b*.  In  this 
case  both  integrals  above  are  close  to  1 as  n -£>oo  . 

2 2 

Finally,  we  note  that  H(b)  is  Schur-convex  in  (b-^,...,bp),  so 
that  further  bounds  on  H may  readily  be  obtained.  In  particular,  it 
can  be  shown  via  such  bounds  that  vn  decreases  exponentially  in  n, 
so  that  we  can  expect  reasonably  good  approximations  for  moderate 
values  of  n. 
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Tab  I e 1_.  Numerical 
for  d = 0 .00(0 .02)0 


Values 


|i(/+id)d;(y) 


20(0.10)4.00  , p= 2(1) 10 ,12 ,15  ,20 


2 

3 

4 

5 

6 

7 

8 

9 

10 

12 

15 

20 

0.00 

. 3333 

. 2500 

. 

. 1657 

. 1429 

. 1111 

. 

.0909 

. 

. 

.0476 

0.02 

. 3418 

.2604 

.2118 

. 1795 

.1567 

.1396 

. 1264 

.1153 

. 1074 

.094$ 

.0814 

.0663 

0.04 

. 3504 

. 2709 

. 2238 

. 1928 

.1710 

.15  4'' 

. 1425 

. I 3 2 • 

.1251 

. 1135 

. 1023 

. 0920 

0.06 

. 3589 

.2815 

. 2 360 

.2063 

. 1857 

. 170  7 

. 1593 

.1506 

. 1437 

. 

. 124$ 

. 1172 

0.08 

. 3676 

.2922 

.2484 

.2202 

. 2008 

.1870 

. 1767 

.1689 

.16  30 

. 1547 

. 1477 

. 1429 

0.10 

. 3762 

. 3030 

.2610 

.2343 

.2163 

.2036 

.1945 

.1877 

.1827 

.1761 

. 1710 

.1681 

0.12 

. 8 6 4 9 

. 3139 

.2737 

. 2486 

.2319 

.2205 

.2125 

.2068 

. 2027 

.1975 

. 1940 

. 1024 

0.14 

. 3936 

. 3249 

. 2865 

. 2630 

. 2478 

.2376 

.2307 

. 2259 

. 2226 

.2168 

.2165 

.21$  7 

0.16 

.4024 

. 3360 

. 2994 

. 2775 

. 26  37 

. 2547 

. 2489 

. 2450 

. 2424 

.2396 

. 2382 

. 2378 

. 

.4111 

. 3470 

. 3124 

. 2921 

. 2797 

.2719 

. 2669 

. 2638 

. 2619 

.2600 

. 2591 

. 2539 

0.20 

. 35  61 

. 3254 

. 3067 

. 2956 

. 2889 

. 2848 

. 2824 

.2810 

. 2797 

.2792 

. 2791 

0.30 

.4636 

. 4136 

. 3901 

. 3786 

. 3730 

. 3704 

. 3692 

. 3687 

. 36  84 

. 3683 

. 3683 

. 3633 

0.40 

.5067 

.4679 

.4524 

. 4462 

.4439 

.4432 

.4429 

. 4428 

.4428 

.442  8 

. 4428 

.4426 

0.50 

. 5487 

.5199 

.5104 

. 5075 

.5068 

. 5066 

. 5066 

. 5065 

. 5065 

.5065 

. 5065 

. 5065 

0 . 6 0 

.5892 

. 5686 

. 56  3 3 

. 56  21 

. 5619 

. 5619 

. 5619 

. 5619 

. 5619 

.5619 

. 5619 

. 5619 

0 . 70 

A i 

.6164 

.6138 

.6134 

.6134 

.6134 

.6134 

.6134 

.6134 

.6134 

,t  : 14 

.6134 

0.80 

.6638 

.6546 

.65  33 

.6532 

.65  32 

. 65  32 

. 6532 

.6532 

.6532 

.6532 

.653: 

0.90 

.6976 

.6918 

.6912 

.6912 

.6912 

.6912 

.6912 

.6912 

.6912 

. 6912 

.6912 

.6912 

1 .00 

. 7288 

. 7253 

. 7251 

. 7251 

. 7251 

. 7 2 1 

. 7251 

. 7251 

. 7251 

. 7251 

. 7251 

. 7251 

1.10 

.7574 

. 7574 

. 7573 

. 7573 

. 7573 

. 7573 

. 7573 

. 7 5 ? 3 

. 7573 

. 7573 

. 7573 

. 7573 

1.20 

. 7836 

.7825 

. 7825 

. 7825 

. 7825 

. 7825 

. 7825 

. 7825 

. 7825 

. 7825 

. 782  5 

. 7825 

1.30  | 

f 

. 8074 

. 8068 

. 8066 

. 8068 

. 8068 

. : 68 

. 8068 

. 8068 

. 8068 

. 8068 

.80$  8 

. 8068 

1.40 

.8290 

. 8287 

. 8287 

. 82  87 

. 8287 

. ■ 87 

. 8287 

.8287 

. 8287 

. =2  87 

. 3287 

1.50 

. 8485 

. 8483 

. 8483 

. 8483 

. 8483 

. 4 

. 8483 

. 8463 

. 8483 

. 8483 

. 8483 

. 8483 

1.60 

. 8661 

. 8660 

. 8660 

. 8660 

. 8660 

. 8660 

. 8660 

. 8660 

.8660 

. 8660 

. 8660 

. 8660 
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